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TECHNICAL PUBLICATION 


LUNAR HABITAT OPTIMIZATION USING GENETIC ALGORITHMS 

1. INTRODUCTION 


The Nation’s Vision for Space Exploration calls for a human return to the Moon by 2020 and a later 
potential human mission to Mars. These future missions are expected to involve much more impressive 
activities than those of the Apollo program. Some missions may last several months, and others may last 
up to 600 days. 1 These long-duration surface missions will require large outposts to accommodate living 
quarters (habitats) as well as indoor laboratory facilities. Transporting the materials required to build the 
necessary habitats will be costly and dangerous. 2 The greatest impediment to extended human presence 
on the Moon or Mars is the threat posed by the harsh environments found on their planetary surfaces and 
on the journey there. The lunar environment is much different than the terrestrial environment; many of 
the load conditions on the Moon are much more severe. 2 

Crew habitats require thermal insulation, structural support, impact resistance, and radiation pro¬ 
tection. The temperature of the habitat and electronics contained within must be controlled. 3 Temperatures 
on the Moon can range from 114 °C at lunar noon to -183 °C during lunar night. 4 The Moon is essentially 
without an atmosphere and Mars’ atmosphere is minimal; therefore, protection against meteoroids is also 
necessary. “Lor astronauts in a pressurized building that is suddenly pierced by a small meteorite, the 
chances of survival are problematic.” 5 Ultraviolet and ionizing radiation penetrates to the lunar surface 
without any attenuation due to the vacuum environment of space and the lack of any lunar atmosphere. 3 
Cosmic radiation and solar particle events pose a risk to the crew on the surface of either the Moon or 
Mars. 


To explore the solar system, a permanent base on the Moon is necessary, and “any large-scale 
permanent colonization of the Moon demands the use of locally available lunar materials, simply because 
the mass of materials needed in construction would be too great and too expensive to import from Earth.” 2 
Through the use of in situ resources and terrestrial materials a habitat can be designed to protect against 
the hazards on the surface of the Moon and Mars. 

The materials chosen for the habitat walls play a direct role in protection against each of the men¬ 
tioned hazards. Multifunctional materials; e.g., radiation shielding that also has good thermal insulation 
properties, are an excellent way to provide protection since materials can be selected that meet not only 
the structural requirements but also afford good radiation and thermal protection as well. 


Choosing the best materials, their configuration, and the amount required is extremely difficult due 
to the immense size of the design space. For example, a three-layer habitat wall created from 20 possible 
materials, all with varying thicknesses, would have the following number of possible solutions: 

(number of layers)*(possible material types)h X)SSlblc thicknesses) _ 72 °° . (1) 

Clearly, some sort of optimization method would be extremely beneficial. 
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2. OPTIMIZATION 


There are several optimization techniques, each having a number of advantages and disadvantages, 
that could apply. Standard optimization techniques are not suitable for problems with the large search 
space that will be encountered when designing a surface habitat. Standard (direct) optimization techniques 
(also referred to as “hill-climbing” techniques) seek local optima by moving towards the direction of the 
local gradient, and indirect techniques search for local extrema by solving a set of (typically nonlinear) 
equations that result from setting the gradient of the objective function equal to zero. Unfortunately, both 
methods only seek local optima. 6 

A genetic algorithm (GA) is a search algorithm based on the principles of evolution and natural 
genetics. A pseudocode example of a GA is shown in figure 1. GAs can be implemented as a population 
of binary strings or as a string of real-valued numbers, each string representing a solution to the prob¬ 
lem. Genetic operators—crossover and mutation—are used to create the succeeding generation from the 
members (strings) of the current population. The algorithm is repeated until a stop criterion is reached; 
e.g., a certain number of generations are processed. 7 Each individual is assigned a fitness based on how 
well it satisfies the objective function, which is calculated through engineering analysis. The probability 
of a design being propagated to the next generation as well as the frequency of occurrence of the design in 
the next generation is directly proportional to the fitness of the individual. 6 


ProcecLreGA{ 

t=0t 

Initialize P©; 

Evaluate P©; 

\Mrile (Not Done) 

{ 

Parents© =Select_Parents(P©); 
Offspring© =ProcreatE(Parents©); 
ExflluatefOffspring©); 

P(t+1)= 

SdectJ5r^vas(P©,Offspring©); 
t=t+t 

} 

} 


Figure 1. Pseudocode version of a GA. 


The use of GAs can greatly enhance the selection of materials that have more than one purpose, 
such as radiation, thermal, and impact protection, and multifunctional materials can save a considerable 
amount of mass when compared to standard design techniques. Combining the multifunctional materi¬ 
als with in situ resource utilization will allow a surface habitat to be designed with a minimal amount of 
up-mass. 
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3. BACKGROUND 


3.1 Objective 

The objective of this work was to develop a cross-disciplinary surface habitat optimization and 
analysis tool that can determine a near-optimal combination of materials for the walls of a surface habitat. 
The tool optimizes the habitat materials for minimum up-mass including the analysis of thermal losses and 
gains, structural integrity, meteoroid impact, and radiation shielding. A pictorial of the various objectives 
is shown in figure 2. 

Thermal Losses Structural Integrity Up-mass Meteoroid Impact Radiation Dose 
Figure 2. Multiobjective parameters for the lunar habitat design tool. 


3.2 Thermal 

The thermal stability of a lunar habitat is of extreme importance to the survival of the crew and 
continued viability of the habitat. A lunar habitat experiences drastic thermal cycling loads during a lunar 
day. Temperature extremes can vary 200 K on the surface of the Moon; therefore, it is necessary to design 
a structure that can sustain such dramatic thermal cycling loads. A “shirt-sleeve” environment is the ideal 
and preferred environment for astronauts. It is common knowledge that active heating and cooling is 
required for lunar habitats, and a reduction of such a need will reduce the up-mass of the required equip¬ 
ment. Therefore, thermal losses are calculated to determine the suitability of a given design. 

3.3 Structural 

Structural integrity of the lunar habitat is necessary to ensure that the habitat does not fail from 
overload conditions. A lunar habitat must be able to withstand many different structural loads including 
residual loads from the habitat itself, the pressure differential of the shirt-sleeve environment and the 
vacuum of space, impacts from meteoroids, and loads during construction phases. All of these loads are to 
be considered when evaluating the structural integrity of the lunar habitat. A safety factor of 1.1 was used 
to analyze a fully loaded habitat. 


3.4 Meteoroid Impact 

A meteoroid is a “naturally occurring solid body, traveling through space, which is too small to be 
called an asteroid or a comet.” 8 It is estimated that meteoroids of =T-mg mass size strike lunar structures 
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yearly. Velocities of meteoroids at the Moon range from 13 to 18 km/s. Meteoroids of =T0 -6 g can generate 
craters of 500 pm in metals. Impacts of larger meteoroids are more rare; however, they pose a much more 
significant threat to the habitat. A meteoroid with a mass of =1 g can produce a centimeter-sized crater in 
metallic materials. 8 


3.5 Radiation 

The radiation environment encountered in interplanetary space consists primarily of energetic ions 
and the galactic cosmic rays (GCRs). Further sources of radiation are solar particle events (SPEs) due to 
solar flare cycles. SPEs are infrequent events that occur with large coronal mass ejections, and result in a 
large amount of particles, mostly protons, moving through the solar system. 9 

GCRs originate outside our solar system, are regulated by the solar wind, and contain high-charge 
ions as well as alpha particles and protons. The heavy ions make up a small portion of the total particle 
fluence (particles per area during a given time); however, they are responsible for approximately half of 
the long-term dose equivalent. 9 

For long-term stays, the crew should have at least 400 g/crn- of regolith covering for adequate 
radiation protection, which does not account for solar flare events. Flowever, 100-200 g/cm 2 of regolith 
intensifies the dosage due to secondary radiation. 10 At least 2.5 m of regolith covering would be required 
to keep the annual dose of radiation at 5 rem. 11 The limit for astronauts in low Earth orbit is 25 rem, and 
the occupational limit for an adult is 5 rem/yr. 

The addition of other materials could decrease the amount of regolith cover required for radiation 
protection. Flydrogen and materials that contain a large amount of hydrogen, such as polyethelyene and 
other polymers, have been shown to be good radiation shielding materials. 9 A combination of lunar rego¬ 
lith and other materials would provide a good barrier to radiation. 

3.6 Materials Database 

The materials database contains all the necessary material properties for each of the analysis 
modules. Currently there are 24 materials in the database, including aluminum, polyethylene, Mylar®, 
regolith cement, polyetheretherketone (PEEK), stainless steel, Kevlar® and glass composites, rubber, and 
ethylene-vinyl alcohol copolymer (EVOH). Other materials can easily be added to the database; the name 
and properties are all that are required. 
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4. COMPUTER CODE DEVELOPMENT 


The computer code was written primarily in MATLAB. Each analysis module was written as a 
separate subroutine. The optimization code is a GA written by the authors at NASA Marshall Space Flight 
Center. 


4.1 Thermal 

4.1.1 Steady State Numerical Thermal Model 

A significant amount of research has been done to develop an analytical model of the multilayer 
insulation heat leak of the lunar habitat under consideration. Analysis of the habitat begins by the applica¬ 
tion of the first law of thermodynamics to the control volume. 

This is a passive system, meaning that no additional heat comes in or is taken away by an active 
thermal control system. A depiction of this system is shown in figure 3, and the corresponding energy bal¬ 
ance equation is 


E m + E a + E, 


out 


( 2 ) 


where 

E- m is the energy into the system from outside sources 

E out is the energy loss from the system 

Eg is the thermal energy generation 

E st is the rate of energy stored within the system. 



E 


out 


Figure 3. Theoretical control volume of the system. 
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The core problem of cooling and/or heating lunar habitats is the thermal radiation due to solar 
heating by direct and indirect means; i.e., directly from the Sun and indirectly from lunar reflection. The 
effective radiation heat flux equation from the Sun (or net energy absorbed from the Sun) is 


GradSun a ' Asun ' '/Sun > 


(3) 


where 

a is absorptivity 

A Sun is the area of the habitat being radiated on by the Sun 
<7Sun ' s the incident heat flux from the Sun. 

The effective radiation from the habitat takes the form of 


QradHab - e ■ O ■ A ■ T s 4 , 


(4) 


where 

e is the emissivity 
a is the Stefan-Boltzmann constant 
A is the surface area of the habitat 
T s is the surface temperature. 

Calculation of the surface temperature of the habitat requires the determination of the roots of a 
fourth degree polynomial. Equating the conduction and radiation heat transfer rates and the polynomial is 
shown in equation (8). 


2cond - Gkad 


AA^=e T ot-«aj - s 4 

^EQ 

(Qtot)^EQ - eaA T s 4 R E Q -T s + 7j =0 
Q ^//eq^Ts 4 0 o -T s (Qjoj)Req +a] 


(5) 
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4.1.2 Finite Element Analysis Thermal Model 

A finite element analysis (FEA) model of the thermal behavior of the lunar habitat was created in 
ANSYS, and an example result is shown in figure 4. 




-0.564x10" 05 165.942 331.884 497.825 663.767 


Figure 4. Example FEA thermal model. 


4.2 Structural 

A structural model was created and written in ANSYS. This structural model includes all load¬ 
ing conditions, including pressure differential and material weight in a 1/6-g environment. An example 
ANSYS plot is shown in figure 5. 



1,709 12,722 23,734 34,747 45,759 56,772 67,784 78,797 89,809 100,822 


Figure 5. Example structural FEA model. 







4.3 Meteoroid Impact 


The meteoroid impact analysis was developed around the Fish-Summers single-plate penetration 
equation: 12 


D v 0.352 „ 0.167 , T , _ ox0.875 
P = Km p p p (Vp cos p) 


(6) 


where 

P = depth of penetration (cm) 

K = material constant 
m p = mass of the projectile 
p p = density of the projectile 
V p = velocity of the projectile 
/3 = angle of impact. 

The projectile density and velocity were assumed to be 0.5 g/cm 3 and 20 km/s, respectively. The 
velocity is conservative, and the density is considered an estimate. The zero-degree angle of impact used 
was considered a conservative assumption. 


The material constant can be approximated for each material with the following equation: 13 


K = 


0.816 

£ 1,1S ^ 


(7) 


where £ = percent elongation and p T = target density. 

The projectile mass was determined to be 0.22 g by iterating between equations (8) and (9) (Rob¬ 
inson, J.H.: Private Communication, April 7, 2005). 14 

r _ ft 0.306 . x^t.38 . , , 2 . 4 x-0.36 ,6 ( . 2x-0.85, ^ 

F — Cq\\C\WI +C 2 ) + C^yWl + C 4 /W + C^lTl ) + C \17l + C^iTl ) } , ( 8 ) 

where 

F = number of particles/m 2 /yr of mass, m, or greater 

c 0 = 3.156 xlO 7 

cj = 2.2xl0 3 

c 2 =15 

c 3 = 1.3xl0- 9 

c 4 = 10 11 

c 5 = 10 27 

c 6 = 1.3 xl0“ 16 

c-j = 10 6 . 
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p _ ~F*A*t 

'NI — e 


(9) 


where 

P NI = probability of no impact 
F = number of particles/m 2 /yr of mass, m, or greater 
A = area (m 2 ) 
t = time (years). 

To be conservative, the mass was calculated assuming a 32-m outer radius, 15-yr time span, and 95 per¬ 
cent probability of no impact. 

Only the outside layer of the habitat wall was analyzed by the impact code because averaging all 
of the material constants ( K) for each layer of the multilayered wall would have resulted in a less accurate 
penetration depth. Furthermore, this forces a tough outer layer to protect the habitat, which is the most 
realistic assumption. 


4.4 Radiation 

The radiation analysis is considered to be a rule-of-thumb, level-zero analysis. The thickness of 
each layer is converted into a corresponding thickness of aluminum and then the dose equivalent is calcu¬ 
lated in the middle of the habitat. 

The 1989 Solar Maximum radiation environment is used in this analysis because it is assumed to 
be a near-worst-case SPE. The radiation model only accounts for protons, but a 5 rem dose is added to 
account for any secondary neutrons created from the primary proton bombardment. 

4.5 Materials Database 

The materials database is contained as a matrix in a separate file accessed by the code. It contains 
all the necessary material properties to run each of the analysis modules. Also, a name matrix is contained 
in the same file, and it contains character strings of the names of each material type. This is used to gener¬ 
ate output files and when creating the analysis plots. 

4.6 Optimization 

A GA was written specifically for the lunar habitat wall optimization problem. This GA uses a 
combination of crossover, mutation, and attrition to search the design region for a lunar habitat optimiza¬ 
tion problem. 

GAs require few control parameters, and assuming the computational time and resources 
are not constrained, settings for these control parameters become less of an issue. However, if high- 
performance GAs that are able to find excellent solutions quickly are desired, then it is important to dis¬ 
cover a set of optimal control parameter settings. In the GA literature, there are basically two approaches 
used in the discovery of high-performance control parameters: (1) An online approach (adaptation) 1?_ 18 
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and (2) an offline approach (discovery of robust control parameters). 19 ’ 20 Online adaptive approaches 
seek to evolve appropriate control parameter settings for the GA during the evolutionary process. Control 
parameter settings that may work well one stage (beginning, middle, or later stages) of genetic search may 
not work well during the other stages. The online adaptive approaches are therefore developed to discover 
the most appropriate parameter settings to apply at the appropriate time during genetic search. 

In the offline approach, the objective is to find a set of control parameter values that work well for 
a set of related applications. Shaffer, et al. and Grefenstette have developed a meta-GA for evolving robust 
control parameters for a simple GA. 19 ’ 20 These parameters have been used as a guide for control param¬ 
eter setting by GA researchers for solving a wide range of problems. 6 In this Technical Publication (TP), 
a meta-GA is used to evolve a set of control parameter settings for population size, mutation rate, and 
mutation amount (perturbation amount or sigma). Each member of the meta-GA was evaluated five times 
to determine a reasonable fitness value for each member of the population. 

4.7 Meta-Genetic Algorithm Results 

The meta-GA focused on searching for three optimal GA parameters: (1) Population size, 
(2) mutation rate, and (3) mutation amount. The plots of these GA parameters are shown in figure 6. The 
optimal GA values determined follow: 

• Population size = 150. 

• Mutation rate = 0.332. 

• Mutation amount = 0.6128. 

The allowable ranges for the GA design parameters were set at 4-200 for population size and 0-1 
for both mutation rate and mutation amount. A plot of the objective function versus generation for an 
example run using the optimal GA parameters is shown in figure 7. 
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Figure 6. GA control parameter plots. 



Figure 7. Meta-GA objective function plot. 
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5. OBJECTIVE FUNCTION FORMULATION 


The habitat optimization tool described in this TP attempts to optimally combine several aspects 
pertinent to a superior lunar habitat design, as follows: 

• Minimize thermal losses. 

• Maintain structural integrity. 

• Provide meteoroid and radiation protection. 

Each of the mentioned parameters is independently calculated based on the material and geometric 
framework of the candidate habitat design. The difficulty lies in the weightings of the individual objectives 
in question. Two things to consider are: (1) Should they all be weighted equally, and (2) how is equality 
measured? These questions must be addressed in order to fully optimize a realistic lunar habitat design. 

The weightings of the objective functions are simply a subjective measure of the designer’s ordered 
priority. It may be the choice of the designer to emphasize the minimization of the active thermal control 
required to maintain the shirt-sleeve environment for the astronauts, or the designer may wish to empha¬ 
size meteoroid and radiation protection. Regardless of the designers choice, it is necessary to scale the 
objectives such that the weightings are initially equal. Recalling that the various objectives are measured 
in heat loss, deflection, penetration depth, and radiation dose. To equally weight these parameters, a series 
of 60 analyses were performed for this study. The objective weightings were changed sequentially to 
equalize the differing objective parameters. The final values are listed in table 1. 


Table 1. Final equally scaled constraint factors. 


Parameter 

Factor 

Up-mass 

Heat loss 

Deflection 
Penetration depth 
Radiation dose 

1 

11.68 

872,227.76 

95,233.43 

51,064.29 


The values in table 1 are simply an equal weighting representation of the individual objectives. At 
this point, the designer has the option of emphasizing one objective over another. Examples of designer- 
weighted objectives are given in section 6. 
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6. RESULTS AND DISCUSSION 


The GA parameters, as chosen by the meta-GA, were applied to a single run of the lunar habitat 
design tool. This final run was compared to the following user-selected GA parameters: 

• Population size = 40. 

• Mutation rate = 0.20. 

• Mutation amount = 0.20. 

In order to provide an equal starting point, the intuitive parameter run utilized the best 40 indi¬ 
viduals from the randomly generated initial population of the optimized parameter run. Figure 8 clearly 
demonstrates that the optimal GA parameters return a better solution with a much quicker convergence. 


5,000,000,000 
4,500,000,000 
4,000,000,000 
3,500,000,000 
3,000,000,000 
2,500,000,000 
2 , 000 , 000,000 
1,500,000,000 
1 , 000 , 000,000 
500,000,000 
0 

0 20 40 60 80 100 120 140 160 

Generation 

Figure 8. Optimal meta- and intuitive-GA objective function comparison. 



6.1 Example Optimization Result 

The results of an example optimization run are shown in figure 9. This configuration was obtained 
with an up-mass weighting of 6; i.e., the calculated up-mass was multiplied by 6 to allow a greater influ¬ 
ence on the final design. Table 2 includes the constraint weightings for this optimization run. The total sum 
of the weightings is 1, not including the up-mass. 
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Habitat Cross-Section Optimal Material Configuration 


Gen: 300 
Q (W): 7,775 
delta (m): 1,0075 
Dose (rem): 12,3318 
P (m): 4,8875 



I-1 Al 1200-H4, t=0.0160 

l l Lunar Cement, t=0.0805 

l l Lunar Cement, t=0.2523 


Table 2. Weightings for example 1 optimization run. 


Constraint 

Weighting 

Up-mass 

6 

Heat loss 

0,3 

Deflection 

0.1 

Penetration depth 

0.3 

Radiation dose 

0.3 


Figure 10 displays the thermal and structural ANSYS analysis plots for final example optimized 

habitat. 
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limum 



r —x 

Maximum \\ 

J u 


89.282 178564 267.817 357.129 

44.611 133.923 223205 312488 40177 


555852 32826 ' 65 , 096 ' 97366 129,635 

16691 48,961 81231 113,501 145,770 


(a) Thermal Analysis (b) Structural Analysis 

Figure 10. ANSYS plots for example 1: (a) Thermal analysis and (b) structural analysis. 


This example demonstrates a design where mostly lunar regolith is chosen, which is expected 
since a large emphasis is placed on minimizing up-mass. 
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7. CONCLUSIONS 


A cross-disciplinary surface habitat optimization and analysis tool has been developed that can 
determine a near-optimal combination of materials for the walls of a surface habitat. The tool optimizes 
the habitat materials for minimum up-mass including the analysis of thermal losses and gains, structural 
integrity, meteoroid impact, and radiation shielding. There is evident value in determining the optimal 
habitat wall configurations. 

The GA tool, as presented in this TP, provides a preliminary optimized habitat to give a design 
team a starting point for more detailed engineering and analysis. The focus of the optimization can be 
changed by varying the weightings of each of the analysis constraints. This adds a great deal of flexibility 
for the user and allows the design to be focused for each application. 

Also presented is a meta-GA for optimizing the GA parameters of the habitat optimization tool. 
This meta-GA has been proven to provide faster results and better solution sets. 
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